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Molecular biology explains function of molecules by their geometrical and electronical structures 
that are mainly determined by utilization of quantum effects in chemistry. However, further quantum 
effects are not thought to play any significant role in the essential processes of life. On the contrary, 
consideration of quantum circuits/protocols and organic molecules as software and hardware of liv- 
ing systems that are co-optimized during evolution, may be useful to overcome the difficulties raised 
by biochemical complexity and to understand the physics of life. In this sense, we review quan- 
tum information-theoretic approaches to the process of DNA replication and propose a new model 
in which 1) molecular recognition of a nucleobase is assumed to trigger an intrabase entanglement 
corresponding to a superposition of different tautomer forms and 2) pairing of complementary nucle- 
obases is described by swapping intrabase entanglements with interbase entanglements. We examine 
possible biochemical realizations of quantum circuits/protocols to be used to obtain intrabase and 
interbase entanglements. We deal with the problem of cellular decoherence by using the theory of 
decoherence-free subspaces and subsystems. Lastly, we discuss feasibility of the computational or 
experimental verification of the model and future research directions. 



I. INTRODUCTION 

According to the central dogma of molecular biology, 
genetic information stored in double-stranded DNA (ds- 
DNA) is duplicated by replication of two strands inde- 
pendently. At each step of the replication, enzyme DNA 
polymerase (DNApol) first recognizes the nucleobase (N 
= {A, T, G, C}) of the template DNA strand. Then, 
it finds complementary of this base (N = {A=T, T=A, 
G=C, C=G}) from the surrounding environment and fa- 
cilitates pairing of these bases through two or three inter- 
base hydrogen bonds. A new dsDNA is synthesized from 
an existing single-stranded DNA (ssDNA) by successive 
pairings of the all nucleobases in this way. 

Conformational changes occurring within DNApol at 
each stage of the replication were demonstrated in de- 
tail by crystallization experiments [38]. Also, all pos- 
sible interactions between amino acid side chains (that 
are likely to be found in the active site of DNApol) and 
unpaired nucleobases were obtained by quantum chemi- 
cal calculations [7]. However, there are still some unclear 
points about the relation of the high fidelity of replication 
process with the base recognition, searching, and pair- 
ing mechanisms [38]. Since active site of DNApo/ that 
contribute to these mechanisms has a particularly com- 
plex structure involving a lot of amino acids [38] . both 
experiments and quantum chemical calculations are still 
insufficient to clarify these mysteries. Thus, until the 
development of more sensitive setups and more sophisti- 
cated calculations, information processing models could 
be useful tools for a better understanding of DNA repli- 
cation. 

During the DNA replication, newly synthesized 
strands elongate with a rate 3,000 nucleotide per minute 
in humans [45] and 30,000 nucleotide per minute in bac- 
teria [32 . Neither DNA binding nor nucleotide bind- 



ing to the DNApol limits this rate, they are very fast 
steps [14] [I5j [38] . Also, replication without proofreading 
and repair mechanisms occurs with an error rate of 10 -4 
to 10 -6 per nucleotide [32] . Such an accuracy can be 
within the constraints of quantum coherent information 
processing. Estimations based on both theoretical mod- 
els [58] and experimental data [TJ [3] give sufficiently long 
decoherence times for the coding nucleobase protons of 
dsDNA [51] . Thus, quantum information processing de- 
scriptions are expected to be explanatory models of DNA 
replication. 

To understand the underlying mechanisms of DNA 
replication several quantum information processing mod- 
els were proposed. For example, Patel [35] formulated 
nucleotide selection from surrounding environment as an 
unsorted database search. He examined the pertinence 
of Grover's algorithm [16 to give an explanation for the 
number of deoxyribonucleotide types used in dsDNA. Al- 
though he successfully modeled base pairing as oracle of 
the algorithm [35] and associated this model with the 
evolution of the triplet genetic code [36 , initiation of 
the search in his model requires the symmetric quan- 
tum superposition of four disparate nucleotides which is 
not quite possible. Wave analogue [37] of this quantum 
search algorithm in which symmetric superposition state 
is replaced by the center-of-mass mode is more realistic 
for enzyme activity. However, if this version of the al- 
gorithm [37] is adopted for the activity of DNApo/, each 
base pairing should begin with the loading of DNApol 
with four different free nucleotides before attempting to 
bind to DNA which is contrary to the present knowledge 

[H [T5l [58]. 

Recently, Cooper [9] [10] modeled base recognition 
mechanism in replication (and transcription) to under- 
stand time-dependent DNA mutations and AT richness 
of DNA. To explain the stability of base pairs, he as- 



2 




sumed that interbase hydrogen bonds are rearranged by 
sequential intermolecular and intramolecular proton tun- 
nelings. In this assumption, interbase tunnelings turn 
bases into their unusual tautomers (N* and in Fig- 
ures [I] [2| pair-by-pair. Then, intramolecular tunnel- 
ings introduce coherent superposition states in which enol 
and imine protons of unusual tautomers are shared be- 
tween two electron lone pairs that belong to a single 
atom. In the recognition step of the model, enzyme 
transcriptase (RNA-dependent DNApo/) makes quantum 
measurements on the coherent states of protons that are 
present on Watson-Crick (WC) edge (Figure [3J. Al- 
though this model is compatible with molecular genetic 
transcription data of bacteriophage T4, some possible re- 
sults of the transcriptase measurements, such as the deco- 
hered states corresponding to tautomers Gq 02 and Gq 00 
[9] [10] do not generate information for any usual tau- 
tomer form: technically, qubit representation of a nucle- 
obase was considered as the tensor product of states cor- 
responding to the presence or absence of each WC edge 
proton. Then, Hilbert space should be 8-dimensional, 
but four bases states, that do not correspond to com- 
mon tautomers, give no meaningful information to the 
enzyme. Therefore, in such situations where the result 
of the measurement is one of these spurious states, en- 
zyme can not recognize the nucleotide base of DNA. This 
expresses an efficiency problem in both recognition and 
searching mechanisms. 



II. ENTANGLEMENT SWAPPING MODEL OF 
REPLICATION 

Estimation of long decoherence times for the WC edge 
protons [31] suggests that nontrivial quantum effects can 




FIG. 2: G-C^G*-C* and G-C-^G*C* tautomeric transitions 
which can be observed [12] [26] [51] [53] in dsDNA. 

be involved in DNA replication. It is easy to show that 
the fastest quantum search algorithm [16 , 57 is only four 
times faster than the slowest classical search algorithm 
for the case of searching complementary nucleotides of 
template bases. Therefore, if quantum coherence is main- 
tained during replication, this should have been evolved 
to increase not only the speed, but also the accuracy. In 
this sense, our model is motivated to investigate quan- 
tum effects increasing both speed and accuracy together 
in the DNA replication. 

A. Intrabase and Interbase Entanglements 

During the DNA replication process free nucleotide 
binds to a solvent exposed pocket within DNApo/ be- 
fore base pairing, whereas template base is flipped out 
of the helix axis and into the active site [38]. Addition- 
ally, it is theoretically known that interaction with water 
molecules can induce transitions to rare tautomer forms 
[T3] . Such higher energetic states can also be mediated 
by interactions with carboxylate and sodium ions [42] 
which are likely to be found both in the solvent exposed 
pocket and active site of the enzyme. Thus, tautomeric 
transitions are likely in both incorporated nucleobase and 
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FIG. 3: Parts of the nucleobases: atoms can be grouped 
according to region they will be found in DNA. Hoogsteen 
(-f-^major groove), Watson-Crick (-H-pairing plane), and Sugar 
(<Hrninor grove) edges are indicated respectively by dashed, 
filled, and plain boxes. Arrows inside the boxes show the 
order of the qubits used in qubit representation. 



template base after recognition and before base pairing. 

In this work, molecular recognition of a nucleobase 
is assumed to trigger a superposition of different tau- 
tomer forms, i.e. \N) W c,i ~> \N)wc,Q = Z)t a *l N *)vrc 
(t = { ,*,#}). Nucleobases A and T have two different 
tautomer forms, whereas G and C have three different 
tautomer forms in allowed transitions (Figures [ll|2j. Su- 
perposition state |N)wc,Q = <^|N)y^c + ^{N/wc re- 
quires the entanglement of first two WC edge atoms of 
nucleobase, whereas |N) wc,Q — ^\~N)wc + <^*|N*)w"C + 
a^|N^)vrc requires the entanglement of all the three WC 
edge atoms. So, superposition of usual and unusual tau- 
tomer forms corresponds to an intrabase entanglement 
of WC edge atoms of the nucleobase. Such interbase en- 
tanglements will be invariant in the situations causing 
tautomeric transitions. Thus, formation of intrabase en- 
tanglements increase the speed of replication if they are 
not fragile to cellular decoherence. 

Possible quantum mechanical transitions AT— ?>A*T*, 
G-C-^-C* and G-C-KT-C* were found in dsDNA by 
density functional theory (DFT) calculations [5QH54] . 
Hence, after base pairing, nucleobase pairs can exist 
in a superposition of states corresponding to different 
tautomer pairs, i.e. |N • N)vfc,o = J2t • N t )wc 

(t = { ,*,#})• However, probability amplitudes for un- 
usual tautomer pairs (/?* and ft) are expected to be 
very small, since transitions to unusual tautomer pairs 
are very rare in dsDNA. Also, nonlocal DFT methods 
[T7HT91 [47] showed that minimum covalent contribution 
to the interbase hydrogen bonds is 38% in AT pairing 
and is 35% in G-C pairing. This can be interpreted as 
quantum mechanical sharing of proton which causes an 



entanglement between the donor and acceptor atoms. In 
this view, superposition states |N • N)wc,o produced by 
the base pairing correspond to interbase entanglements. 

We propose that recognition triggers two two-particle 
entanglements in the \N)wc,Q states of A and T and two 
three-particle entanglements in the |N) wc,Q states of G 
and C (Figure [4]). If we consider each interbase hydro- 
gen bond as a two-particle entanglement, there are two 
two-particle entanglements in | A-T) wc,o state and three 
two-particle entanglements in \G-C) wc,o state after base 
pairing. In order to turn intrabase entanglements into in- 
terbase entanglements, recognition should be followed by 
an irreversible transformation. Therefore, in our model, 
base pairing is described by a multiparticle entanglement 
swapping [4] in which DNApol swaps intrabase entangle- 
ments with interbase entanglements (Figure [4]). 




FIG. 4: Entanglement swapping model of replication. Each 
circle represents the WC edge of related nucleobase. Atoms 
on WC edge are shown by small dark points. A line link- 
ing two of such points means that there is an entangle- 
ment between the atoms shown by these points. Recog- 
nition process is described by a unitary transformation U: 
\N)wc,i — > \N)wc,Q, whereas base pairing is described by an 
irreversible transformation S: |N) wc,q — > \N)wc,o 



While van der Waals interactions between nucleotide 
bases in ssDNA were modeled by entanglement in a re- 
cent study [40] , present study is the first example of mod- 
eling hydrogen bonds by entanglement and using (multi- 
particle) entanglement swapping in living systems. Intu- 
itively, biomolecules appear to be classical objects since 
their de Broglie wavelengths are comparatively smaller 
than their actual size due to their huge molecular mass 
and high temperature. However, it is both theoretically 
and experimentally shown that entanglement can occur 
in macroscopic and hot non-equilibrium systems, such as 
biological ones [2] H [48] [49] . 
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B. Qubit Representation of Input and Output 
Nucleotide States in Replication 

Recognition process requires formation of at least two 
hydrogen bonds between amino acid side-chains of the 
DNApol and nucleobase [TJ [43]. Such pairs of hydro- 
gen bonds can occur over one of the three parts of nu- 
cleotides [7] shown in Figure [3j In consensus, hydrogen 
bond donor and acceptor atoms of bases are only O and 
N atoms. However, there are a small number of computa- 
tional observations in which C atoms of nucleotide bases 
have the ability to make blue-shifting hydrogen bonds 
[27] . In this respect, when electronic configurations of 
the individual O, N, and C atoms on Hoogsteen (if), 
Watson-Crick (WC), and Sugar (S) edges (Figure |3| are 
considered, it is found that each atom has two different 
energy states: a relatively lower energy state as acceptor 
and a relatively higher energy state as donor (Figure [5]). 
These lower and higher energy states can be regarded as 
qubits |0) and |1), respectively (Figure |5j. Then, reli- 
able qubit representations can be written down for all 
the three edges of each nucleobase (Table [I]). 



TABLE I: Qubit representations of usual and unusual tau- 
tomer forms found in the allowed transitions (Figures [I] [2]): 
|0) and |1) states are assigned according to the absence and 
presence of a proton that can be shared in a hydrogen bond 
and order of the qubits are determined as shown in Figure [3] 
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C. Quantum Aspect of the Enzyme Action 

If states of the nucleobases which are measured by 
DNApol live in a Hilbert Space whose dimension is larger 
than the number of these states, there can be unavoid- 
able efficiency problems in both recognition and search- 
ing mechanisms. Also, states corresponding to usual tau- 
tomer forms should be orthogonal to each other. Under 
these conditions, DNApol should recognize bases in both 
free nucleotide and ssDNA cases only over the H edge 
according to Table [I] It is known that sequence-specific 
dsDNA binding proteins usually interact with the major 
groove [41] and so, they recognize nucleobases of dsDNA 
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FIG. 5: Electronic configurations and qubit representations 
of the O, N, and C atoms: higher energy state |1) (lower 
energy state |0)) corresponds to the presence (absence) of a 
proton which is bonded to that atom and which participates 
in interbase hydrogen bonds (Figures [l] [5]). Configuration 
indicated by * is not present in any tautomer form. However, 
it is possible to observe it in blue-shifting hydrogen bonds of 
DNA. 



over the H edge, too. Such a coincidence is not a sur- 
prise, since reading information from dsDNA and ssDNA 
by different proteins ought to be based on similar princi- 
ples. 

If DNApo/ makes a quantum measurement on the state 
|N)h to recognize a nucleobase, first qubit gives infor- 
mation about purine-pyrimidine distinction, whereas the 
second one gives information about imino-enol distinc- 
tion. In this sense, DNApo/ should pair bases whose 
qubit representations are complementary to each other 
(see Table |l|. Not only correct base pairings, but also 
mispairings like AC* and G*T can be accounted for by 
this assumption. 

A quantum measurement requires an entanglement be- 
tween the measuring device and measured system. In this 
case, it can be considered as a hydrogen bonding between 
the DNApo/ and the nucleobase. Also, it is reasonable to 
assume that entanglement between the DNApol and the 
nucleobase should be maximal since an accurate measure- 
ment requires strong coupling between measuring device 
and measured system. 

Hypothetically, a proton transfer between the DNApol 
and the second atom of ii" edge (or equally the first atom 
of WC edge) which occurs during the recognition, can 
trigger a tautomeric transition (Figure [6J. Since such a 
transfer has a quantum nature in a hydrogen bonding, 
recognition can trigger a transition to the superposition 
of usual and unusual tautomer forms by a unitary trans- 
formation U. This mechanism is a toy model of evolution 
of the basis states \N)wc,i (Table [T]) into superposition 
states \N) WC ,Q = Et at \^ )wc,i (t = { ,*,«}) as fol- 
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FIG. 6: An hypothetical mechanism for the tautomeric tran- 
sitions of nucleobase G by proton transfer between enzyme 
DNApo/ and nucleotide. Grey structure represents the active 
site of DNApol and arrows show proton and electron dereal- 
izations. 



formation of intrabase entanglement in N breaks the 
maximal entanglement between DNApo/ and N. Then, a 
maximal entanglement between DNApol and N becomes 
possible. 

Similarly, recognition of N induces an intrabase entan- 
glement in N which disentangles DNApoZ from N. This 
disentanglement allows DNApo/ to bond N and N to- 
gether and then to bind to the subsequent N of ssDNA 
in a quantum mechanical way. Therefore, formation of 
intrabase entanglements not only prevents the uncontrol- 
lable tautomeric transitions caused by cellular environ- 
ment, but also provides separation of DNApol from one 
nucleotide and binding of it to another. 

After base pairing, nucleobase pairs should exist in a 
superposition of states corresponding to different tau- 
tomer pairs. Since state of a hydrogen bonded atom pair 
can be written as the Bell state \J3 01 ) = (|01) + \10))/y/2, 
these superposition states |Ni • N 2 )v^c,o are actually in- 
terbase entanglements. Therefore, it can be said that in 
the case of G-C pair, there are three two-qubit entangle- 
ments in |/?oi) state, and in the case of AT pair, there 
are two two-qubit entanglements in |/3 i) state. In order 
to turn intrabase entanglements into interbase entangle- 
ments, U should be followed by an irreversible transfor- 
mation S. Thus, in our model, base pairing occurs as 
a multiparticle entanglement swapping in which DNApo/ 
swaps intrabase entanglements with interbase entangle- 
ments (Figure 



D. Quantum Circuit for Intrabase Entanglement 
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Now, it is clearly seen that these superpositions of dif- 
ferent tautomer forms are nothing else than intrabase 
entanglements of the atoms on WC edge. Nucleobases 
A and T have two different tautomer forms, whereas G 
and C have three different tautomer forms in allowed 
transitions (Figures [TJ |2j. Thus, after recognition, we 
observe two two-qubit entanglements in the states of A 
and T, while there are two three-qubit entanglements in 
the states of G and C (Figure [4]). 

Recognition process of complementary nucleobase N 
also involves a quantum measurement in which a max- 
imal entanglement is formed between DNApo/ and N. 
However, DNApol can not bind to N in such a quantum 
mechanical way until it disentangles itself from N. This 
is because of the entanglement monogamy (or polygamy) 
[6l [24] which roughly says that if A and B are maxi- 
mally entangled, then any one of them can not be simul- 
taneously entangled with C. In the context of monogamy, 



a) 



SP(6) 



sp(s 



-e- 



b) 



SP(6) — SP'(O) 



43 sp(0) 



-©- 



O- 



FIG. 7: Two of the possible quantum circuits for transforma- 
tion U which turns \N)wc,i states into the \N)wc,q states. 
Superposition matrices SP(9) and SP'(0) of controlled — 
Superposition gates equals to the multiplication of rota- 
tion matrix R(6) and Pauli-Z matrix with different orders: 
SP(0) = R(6) x Z and SP' (0) = Z x R(0). 

Two candidates for the transformation U are shown in 
the Figure [7J To provide a decoherence-free (DF) sub- 



6 



system [2TJ |23l [29j [30] , we take the angles and <p in the 
second quantum circuit (Figure |7]-b) as arccos(V r L/v // 3) 
and arccos(l/v / 2), respectively. Then, \N)wc,Q states 
are obtained as: 
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To consider each base pair as an intact system, tensor 
products of these states should be taken. 
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We reorder qubits of these product states in such a 
way that hydrogen bonded atom pairs come next to each 
other in order to clarify base pairing. Then, we get 
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E. Swapping Protocol for Interbase Entanglement 

Swapping intrabase entanglements to interbase entan- 
glements can be achieved by a three-step protocol S as 
follows: 

1. Reordered base pair states are subjected to a trans- 
formation V as shown in Figure [12) Then, fifth and sixth 
qubits of the G-C (or C-G) pair become |01), whereas fifth 
and sixth qubits of the AT (or T-A) pair become |11). 
After this transformation, any improper base pair exists 
in a superposition of states in which third qubit pair is 
always |00) or |10). 

2. If the sixth qubit is |0), first and second qubit 
pairs undergo a transformation A as shown in Figure 
[9]-a. Otherwise, these qubit pairs are transformed with 
transformation B (Figure [9}b). Then, first and sec- 
ond qubit pairs of proper base pairs collapse into Bell 
state |/3oi) = (| 01) + 1 10)) /a/2, whereas first and second 



FIG. 8: Quantum circuit for S which swaps intrabase entan- 
glements to interbase entanglements: H is the Hadamard 
gate which equals to SP(tt/4) denned in Figure [7] Transfor- 
mations A and B are shown in the subsequent figure. See 
Appendix for the details of transformation V. 



qubit pairs of improper base pairs collapse into Bell state 

|/3 11 ) = (|01)-|10))/v / 2. 

3. Firstly, fifth qubit is passed through a Hadamard 
(H) gate. Then, sixth qubit is converted by NOT (X) 
gate if the fifth qubit is |1). After this step, third 
qubit pair of the G-C (or C-G) pair becomes |/3oi)> 
whereas third qubit pair of the A-T (or T-A) pair be- 
comes |/3n). In contrast, fifth and sixth qubits of any 
improper base pair exists in one of the Bell states |/?oo) = 
(|00) + |ll))A/2 or |ft > = (|00> - |ll))A/2. 
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FIG. 9: Quantum transformations A and B used in S: they 
are actually modified Bell measurements. Modifications made 
by X which is the NOT gate (or Pauli-X gate) and X's super- 
script (Mi or M2) is the outcome of the measurement which 
is done immediately before it. a) A transforms any Bell state 
\f3ij) into the Bell state \/3n). b) B transforms any Bell state 
into the Bell state |/3oi)- 

Immediately after S, proper |Ni • N2 }wc,o states are 
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written in terms of the Bell states as follows. 

|A • TVc, =|T • A) wc ,o = |/3oi)|A>i)|/3ii), (5) 
|G • C) wc ,o=\C ■ G) wc ,o = |A)i>|A)i)|A)i>. 



F. Biochemical Realizations of Quantum 
Circuits /Protocols 

Pauli-X transformation of controlled- NOT gates used 
in the quantum circuit of U (Figure [7]-b) converts 1 1) at 
into \0)n- When state of the DNApo/ is also consid- 
ered with the subscript E, this transformation should be 
\1)n\Q)e \Q}n\1}e' Since |0) and |1) states of an atom 
respectively correspond to the absence and presence of a 
proton bonded to that atom, this transformation can be 
regarded as a proton tunneling from the nucleobase to 
DNApol through the atom on which the gate acts. Vice 
versa is possible for the action of Pauli-X transformation 
on the state \0}n- 

Other gates used in the quantum circuit of U are 
controlled- SP and -SF gates. When argument of SP 
transformation equals to arccos(l/\/2), it transforms 
\0) N into (|0) + |l))/\/2 and \1) N into (|0) - |l»A/2. The 
former transformation should be \0)n\1)e — >• \Poi) ne by 
taking into account also the state of DNApo/, whereas the 
latter transformation should be \1)n\0)e —> \Pu)ne- So, 
the action of the SP transformation on the state \0)n can 
be considered as formation of a hydrogen bond between 
the nucleobases and enzyme through the atom on which 
gate acts. On the contrary, an antibonding should occur 
by the action of SP transformation on the state |1)at. 
This is because of the fact that free energy in the state 
\Pn) ne is greater than the one in which there is no in- 
teraction. Since entanglement measure of the generated 
state changes when the argument is changed, action of 
SP transformation can produce bondings/antibondings 
with different strengths for different arguments. Action 
of SF transformation on \0)n produces the same state 
as the SP action on \1)n an d vice versa. Thus, action 
of SF transformation can be also considered as bond- 
ings / antibondings. 

Both the proton transfer and hydrogen bonding are 
the usual tasks done by enzymes and there are some ev- 
idences for the unignorable role of quantum effects and 
dynamics on the enzymatic reactions [25 , 28 , 44 . There- 
fore producing an intrabase entanglement by transfor- 
mation U is a possible action performed by the enzyme 
DNApol (Figure [10 ). 

Besides U, swapping protocol S also includes 
controlled- NOT, -SP, and -SF gates which are regarded 
as respectively proton tunneling and hydrogen bond- 
ing/antibonding between a nucleobase and DNApol. Ad- 
ditionally, S contains swap gates (Figure [l2|) and modi- 
fied Bell measurements (A and B). Swap gate exchanges 
the states of two qubits on which it acts: 1 10) — » 1 01) and 
1 01) — >• 1 10). So, it can be interpreted as a proton tun- 
neling similar to interpretation of Pauli-X transforma- 
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FIG. 10: Simple depictions of the actions of the NOT, SP, 
and SF gates: atoms whose energy states are represented by 
qubits are the ones inside the shaded spheres. Evolution of 
the \N)wc state of the nucleobase due to the action of gate is 
shown from left to right and |+) qubit equals to (|0) + 1 1))/ y/2 
if equals to arccos(l/v / 2)- 
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FIG. 11: Simple depictions of the actions of the swap gate and 
modified Bell measurement B: atoms whose energy states are 
represented by qubits are the ones inside the shaded spheres. 
Evolution of the reordered |N • ~N)wc state of the nucleobase 
pair due to the action of gate is shown from left to right and 
|An) Bell state equals to (|01) + |10))/^2. 



tion of controlled- NOT gate. However, atoms on which 
swap gate acts can belong either to the same nucleobase 
or to the different nucleobases in the base pair. Hence, 
this proton tunneling should be considered inside a nu- 
cleobase or between the two different bases. 

Outcome of the modified Bell measurement B is the 
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Bell state |/?oi)> whereas the outcome of the modified 
Bell measurement A is the Bell state |/3n). Thus, these 
measurements can also be thought as a hydrogen bond- 
ing/antibonding. Contrary to the ones produced by SP 
and SF gates, this bonding/antibonding is between the 
two nucleobases and its strength is always maximum. 
Consequently, S consists of nothing more than proton 
tunneling and hydrogen bondings/antibondings which 
are the usual tasks done by enzymes like DNApo/ (Fig- 
ures 



10 11). 



Immediately after the entanglement swapping, states 
of proper base pairs are found as in Equation [5] Ac- 
cording to these states, AT (or T-A) base pair has two 
hydrogen bonds and GC (or C-G) base pair has three 
hydrogen bonds as in the actual case. However, these hy- 
drogen bonds have a maximum strength since Bell states 
are maximally entangled. Amplitudes in these states 
should change by the quantum evolution in the presence 
of the asymmetric double well potentials of the hydro- 
gen bonded atom pairs. Thus, strengths of the hydro- 
gen bonds should gradually decrease to the actual ones. 
Moreover, there is an antibonding between the last atom 
pair in A-T (or T-A). These atoms repulse each other 
because of the higher free energy of antibonding, but 
strength of this repulsion should also decrease by time. 
Since one of the atoms in this antibonding is C atom, 
final strength of the repulsion should be negligible. 

On the other hand, both first and second atom pairs 
of the improper base pairs have an antibonding after the 
entanglement swapping. Final strength of the repulsions 
due to these antibonding interactions are not negligible 
and so, they should destabilize and separate the improper 
base pairs. However, state of the last atom pair in these 
base pairs are obtained as |/3qo) or |/?io). Since total 
proton number of the base pair does not remain constant 
after collapsing to these states, atom pair and DNApol 
can not separate from each other. In fact, these Bell 
states should be treated as an entanglement between the 
atom pair and DNApo/ when state of the enzyme is also 



under consideration: 100)^^^ (ll)^^ 



-Hoo) 



NN 



11)r± 



|ll)jyjy|00) e- It can be that it is this entanglement which 
keeps DNApo/ in place till the correct N comes along. 
Both of the asymmetric potentials and destabilization of 
the base pair should weaken this entanglement. When 
entanglement is weakened enough, DNApol can bind to 
the correct N because of the converse monogamy [20] 
which roughly says that if A and B are weakly entangled, 
then any one of them could be strongly entangled with 
C. After that, a Pauli-X transformations can fix the total 
number of protons on the improper base pair and make 
incorrect N separable from the complex. 

Neither U nor S is unique for the given model. How- 
ever, this is not a disadvantage since there are several 
DNApo/ species and families with different replication fi- 
delities. This diversity in replication fidelity of DNApol 
can be accomplished by different U and S pairs. 



G. Effects of Cellular Decoherence 

The intact system which is exposed to decoherence is 
the whole nucleobase - DNApol complex. Hence, states 
of the nucleobases alone are not sufficient to determine if 
decoherence has a significant effect on the transformation 
U or does not. To draw a complete picture of interac- 
tion, assume that there are q hydrogen bond acceptors 
and (k — q) hydrogen bond donors in the active site of 
DNApo/. If so, enzyme's active site can be represented 
by the state \0)® q ^ after a proper ordering 

in which all |0) qubits are put to left of all |1) qubits. 
Then, we can obtain the initial state of the nucleobase - 

DNApol complex as = \N) wc ,i <8> \0)% q <8> \l)f h ~ q) . 
Cellular decoherence effect on the state \s) can be simpli- 
fied as a weak collective decoherence (2TJ [30] which turns 
|1) states into e^|l), while |0) states remain unchanged. 
Since we have already considered |0) and |1) states of an 
atom respectively as the absence and presence of a pro- 
ton bonded to that atom, this simplification makes sense: 
decoherence can not affect an absent proton. 

When spacing between the qubits is smaller than the 
wavelength of the radiation field which acts as a boson 
bath for the qubit system, collective decoherence domi- 
nates among the others [2lJ [29] [30] . In this sense, weak 
collective decoherence seems to be relevant not only to 
long-range electrostatic interactions with the intracellu- 
lar ions, but also to short-living couplings with the ther- 
mal reservoir. Investigation of the effect of thermal reser- 
voir on the internal DNA mobility requires a lattice dy- 
namic approach based on an atomistic description of the 
molecule [56]. According to the appropriate methods 
given in [22, 55 , maximum frequency of the vibrational 
modes in DNA is a few hundreds of cm -1 at room tem- 
perature. Corresponding phonon wavelength is in the or- 
der of fim and this is quite longer than the qubit spacing 
in our model, which is no more than 3 A. Thus, phonon 
bath can not distinguish the qubits and collective de- 
coherence is expected to be the dominant decoherence 
mechanism in the DNA replication. 

Defining a variable \k which equals to the number of 
|0) qubits minus the number of 1 1) qubits in a state over 
K qubits, [2TJ [30] showed that subspaces of Hilbert space 
spanned by the states with constant \k are decoherence- 
free (DF) during a collective dephasing process as de- 
scribed above. Then, a DF subspace for a specific A^ is 
denoted as DFS k (Xk) EDGES- 

In our model, there is no proton exchange between the 
system and its environment during the transformation U. 
For example, when a |1) qubit of nucleobase turns into |0) 
qubit after a controlled- NOT gate, a |0) qubit of DNApo/ 
should also turn into |1) qubit, since Pauli-X transforma- 
tion of cont r oiled- NOT gate corresponds to a proton tun- 
neling between the nucleobases and DNApo/. Therefore, 
value of the \K=3+k remains fixed and transformation U 
does not take any state \s) out of the DFSz+k(2q — fc + 1) 
or DFS3+k(2q — k — 1). This means that decoherence is 
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avoided during the formation of intr abase entanglements. 

Each state that corresponds to an intrabase entangle- 
ment (Equations [I] [2| lives in one of the DF subspaces 
DFS 3 (+1) = Span{j001), |010), |100)} and DFS 3 (-1) = 
Span{|011), 1 101) , 1 110) } . However, intrabase entangle- 
ments should immune not only to the weak collective 
decoherence, but also to the strong collective decoher- 
ence: formation of the intrabase entanglement separates 
the recognized nucleobase from the DNApol This sepa- 
ration removes the isolation provided by the enzyme and 
exposes the nucleobase to the cellular environment. As- 
sumption of the collective decoherence to be weak may 
lose its validity by the removal of the isolation. Thus, ef- 
fects of the strong collective decoherence on the intrabase 
entanglements should also be investigated. These effects 
can be understood in terms of the actions of Pauli-X, 
-Y, and -Z transformations [2TJ|30]. 64 x 64 Pauli spin 
matrices transform intrabase entanglements as follows: 

S x : \N)wc,q\N)wc,q -> \N)wc,q\N)wc,q (6) 
S y : |N)^c,q|N)^c,q -> \N)wc,q\N)wc,q 
S z : \N)wc,q\N)wc,q -> \N)wc,q\N)wc,q 

Since transformation S produces same interbase entan- 
glements from |N)^c,q|N)w^*,q and \N)wc,Q l N ) wc,q 
states, effects of the strong collective decoherence on the 
intrabase entanglements seem to be trivial. This allows 
DNApoZ to safely search for the complementary free nu- 
cleobase after the recognition of nucleobase of ssDNA 
and to safely continue pairing of bases after finishing the 
search. 

Conservation of the proton number of nucleobase - 
DNApol complex is trivial under the actions of swap 
gates and modified Bell measurements. Hence, decoher- 
ence suppression during and after the swapping protocol 
S can be demonstrated in a similar way as is done for the 
transformation U. 



III. DISCUSSIONS 

In the presence of strong collective decoherence, small- 
est DFS (DF subspace or subsystem) in which at least one 
qubit of information can be encoded is a three qubit sub- 
system DFS K = 3 (J = 1/2) [2T1I30] . Construction of this 
subsystem involves the use of four three-particle J = 1/2 
states. So, physical implementation of the computation 
inside the smallest DFS requires the use of four distinct 
building blocks each of which participates in the compu- 
tation over three physical qubits. Number of orthogonal 
states living inside the smallest DFS imposes a restriction 
on the number of different building block types, whereas 
required qubit number for the computation inside the 
smallest DFS puts a limit on the atom number present 
in the interaction region of these building blocks. 

Two of the four three-particle J = 1/2 states that are 
constructed by adding a two-particle J = state to a 
one-particle J = 1/2 state imply the entanglements of 



first two physical qubits of the corresponding building 
blocks. Likewise, other two three-particle J = 1/2 states 
that are constructed by adding a two-particle J = 1 state 
to a one-particle J = 1/2 state imply the entanglements 
of all the three physical qubits of the remaining build- 
ing blocks. Also, two states constructed by the same 
way should be associated in an appropriate way to en- 
code a DF qubit. The conventional way of DF encoding 
inside the subsystem DFSk= 3 (J — 1/2) is taking the 
superpositions of these states. However, quantum super- 
position of the corresponding building blocks are unlikely 
to be formed. Pairing of the building blocks, whose en- 
tangled qubit numbers coincide, would be a more rea- 
sonable way in a biological sense. Such a pairing can 
be obtained by swapping intramolecular entanglements 
with intermolecular entanglements. Since first two of the 
four J = 1/2 states exclude the last physical qubits from 
the intramolecular entanglements, intermolecular entan- 
glements obtained by the pairing of the corresponding 
building blocks should not include the entanglement of 
their last physical qubits. On the contrary, pairing of 
the other two building blocks should result in three in- 
termolecular entanglements between them. 

In fact, if single-particle \j = 1/2, = 1/2) state 
is represented by |0) qubit, four three-particle J = 1/2 
states are nothing else than the \N)wc,Q states corre- 
sponding to intrabase entanglements. It would be in- 
teresting if we could bring the similarities between the 
discussed scenario and processing of genetic information 
beyond an analogy. 

Although tens of nucleotide derivatives are found in na- 
ture, especially in tRNAs, genetic information is encoded 
by only four nucleotides in the dsDNA. These deoxyri- 
bonucleotides usually form Watson- Crick base pairs in 
which pairing occurs over the three WC edge atoms of 
the nucleobases (Figure [3|: A and T form a pair through 
two interbase hydrogen bonds, whereas G and C form a 
pair through three interbase hydrogen bonds. The third 
atom in WC edge of adenine is a carbon atom which is 
responsible for the lack of one hydrogen bond in the AT 
pair. Preference of the carbon rather than a more elec- 
tronegative atom as the third WC edge atom of adenine 
may be a coincidental event stand out in the evolution of 
nucleobases. Such a restriction on both nucleobase types 
and interbase interactions used in the usual structure of 
the genetic information should have an explanation in 
the evolutionary basis. 

One of the essential problems of any information pro- 
cessing in living systems is the presence of unavoidable 
noise. To survive, organisms must have some special 
computational strategies for coping with this problem. 
Since the nature of molecular realm is quantum mechan- 
ical, making computation inside noiseless DFS may be a 
favorable strategy. However, organisms must process in- 
formation not only more accurately, but also more pow- 
erfully. So, any discovery of biological DF computation 
should have been followed by further optimizations on the 
use of resources like energy, time, and number of physical 
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building blocks. Such a resource optimization actually 
corresponds to a minimization of the qubit number re- 
quired in the computation. In this sense, restriction on 
both nucleobase types and interbase interactions used in 
the usual structure of the genetic code may have been a 
result of the further resource optimization in biological 
DF computation. 



IV. CONCLUSIONS 

Since all of the steps in both U and S can be expressed 
as proton transfer and hydrogen bonding, the replication 
scenario proposed in our model could be tested step by 
step with the help of computational methods of quantum 
chemistry. 

In addition to computational tests, some experimental 
setups may be designed to explore some of the predic- 
tions of the model. For example, states of base pairs im- 
mediately after the pairing are obtained as in Equation 
^1 Indeed, |Ni • N2 }wc,o states are superpositions of dif- 
ferent tautomer forms with equal probability amplitudes. 
As discussed before, these amplitudes should change by 
the quantum evolution in the presence of the asymmet- 
ric double well potentials of the hydrogen bonded atom 
pairs until they reach the actual values. Then, it can 
be hypothesized that if this evolution can be prevented 
in a proper way, probability of point mutations due to 
the formation of rare tautomer forms should be higher 
than the ones obtained in vivo. This may be achieved 
by sufficiently decreasing the time periods between two 
successive replications. 

However, there are some deficiencies in the description 
of the present model which should be removed before any 
computational or experimental test. Most important de- 
ficiency is the absence of enzyme's state in the compu- 
tations. To obtain a more realistic description, a state 
should be assigned to the enzyme and the whole process 
including U and S should leave this state invariant at 
the end. In |1991[ a self-consistency condition for a quan- 
tum state was introduced to describe and understand a 
disparate interaction [11 . A similar utilization of the 
self-consistent states in the description of the enzymes 
seems to be appropriate. This is possible, even though 
the state of an enzyme is expected to be mixed, and that 
the purification of a mixed state in |Deutsch| s formalism 
[TT] is impossible [39 j. We note that proofs given in [39] 
are not valid for the mixed states which are used to de- 
scribe enzymes under the self-consistency condition. We 
plan to investigate this further in the near future. 

Finally, entanglement swapping may be a basic tool 
used by enzymes and proteins in the cellular environ- 
ment. If so, similar models may be developed for amino 
acid - tRNA, aminoacyl-tRNA - mRNA, and amino acid 
- amino acid interactions in the protein synthesis. If suc- 
cessful models for these interactions can be developed, 
then we can achieve a deeper understanding of the role 



of the quantum effects and dynamics on the cellular in- 
formation processing. Perhaps, entanglement swapping 
will join and contribute to the debates on the universal 
triplet genetic code [35j [36j [46] and on the mechanism 
behind adaptive mutation (3TJ [33l [34] after such models. 

All in all, quantum effects are used mainly for the 
determination of molecular shapes, sizes and chemical 
affinities in molecular biology and biochemistry. Al- 
though functions of bio-molecules are explained by struc- 
ture, such as the complementary geometries of molecules 
and weak intermolecular hydrogen bonds in nucleobase 
pairs, further quantum effects are not thought to play 
any significant role in the present biochemical complex- 
ity. However, they may be more useful tools to under- 
stand the physics of life if quantum circuits/protocols and 
organic molecules are considered as software and hard- 
ware of the living systems. Reconsideration of evolution 
as co-optimization of hardware (structure) and software 
(function), reconciles two opposite approaches: natural 
selection and self-organization. Thus, emergence of the 
life as a biochemical complexity may be demystified in 
the context of quantum information theory. 
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Appendix A. Transformation V used in S 

Transformation V used in S has thirty seven gates as 
shown in Figure [l2j Such a number seems to be very 
large for an efficient replication process. However, Fig- 
ure [12] is a general representation and all of the gates 
are not effectively used in each base pair (see Table [il]). 
In fact, average effective gate number is approximately 
seventeen. 



TABLE II: Effective gate number of transformation V (Figure 
12) for each base pair Ni- N2. 
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FIG. 12: Quantum transformation V used in S: quantum gates indicated by a double x are swap gates, a) First part of the 
transformation V which makes the states of the proper base pairs orthogonal to each other and to the states of the improper 
base pairs, b) Second part of the transformation V which converts the third qubit pair into 1 01) in G-C (or C-G) pair and into 
1 11) in AT (or T-A) pair. In the case of an improper base pairing, it produces a superposition of states in which third qubit 
pair is always 1 00) or |10). 
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